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ORTHOGONAL POLYNOMIALS OF CONTINUOUS AND 
DISCRETE VARIABLES 
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Abstract. We develop general expressions for the raising and lowering operators that belong 
to the orthogonal polynomials of hypergeometric type with discrete and continuous variable. We 
construct the creation and annihilation operators that correspond to the normalized polynomials and 
■ ,■ ■ study their algebraic properties in the case of the Kravchuk/Hermite Meixner/Laguerre polynomials. 

'NT . 

o ■ 

■ 1. Introduction. In a previous paper [1] we have developed a method to con- 
I struct raising and lowering operators for the Kravchuk polynomials of a discrete vari- 

f— I ■ able, using the properties of Wigner functions, and to calculate the continuous limit 

I to the creation and annihilation operators for the solutions of the quantum harmonic 

• oscillator. 

^sO I In this contribution we apply the same method to other orthogonal polynomi- 

als of discrete and continuous variable. We give general formulas for all orthogonal 
" ^ I polynomials of hypergeometric type [2]: difference/differential equations, recurrence 

■ relations, raising and lowering operators. 

I With the help of standard values we calculate these equations for the normalized 

. functions of Kravchuk- Wigner and Meixner-Laguerre polynomials, and we construct 

' the corresponding creation and annihilation operators. 

. The motivation of this work is to implement the study of the Sturm-Liouville 

I problem in continuous case with the discrete one, in particular the connexion between 

the eigenfunctions and the creation and annihilation operators [3] and [4]. This ap- 
proach is becoming very powerful in the lattice formulation of field theories, where 
the physical properties of the model are analyzed in the lattice before the continuous 
I limit is taken [5] 

I 2. Basic relations between orthogonal polynomials of continuous and 

discrete variable. A polynomial of hipergeometric type of continuous variable sat- 
^ I isfies the following fundamentals equations 

■'^ ' 1. Differential equation: 

^ : (Cl) a (s) y\ (s) + T (s) y'„ (s) + A„ y„ (s) = 

where a (s) and r (s) are polynomials of at most second and first degree, respectively 
and Xn is a constant. This differential equation can be written in the form of an 
eigenvalue equation 

(cr (s) p (s) ?/'„ (s))' + A„p (s) 2/„ (s) 

where p{s) is the weight function satisfying (f7(s)p(s))' — t (s) p (s) and A„ = 
-n(r' + ii^CT"). 

2. Orthogonality relations: 

b 

Un (s) Vm (s) P (s) ds = Snmdl 

with dn some normalization constant, 
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3. Three term recurrence relations: 

(C2) sy„ (s) = anVn+i [s) + (inVn (s) + Tn^/n— 1 



where a„, /3„, 7„ axe constants. 

4. Raising and lowering operators: 



(C3) 
where 



<^ (s) Vn (s) = 



An 



Tn [s) Vn (s) 



Bn 



+ 1 



r„ (s) = T (s) + na' (s) 

A2n+1 



r'n = t' + na" 



2n + l 



We can modify formula (C3) in a more suitable form. 
Prom 



n-l . ^ \ 

a„ = B„ ^r' + -(n + fc-l)a"J , ao = Bq 



k = 

we can prove the following identity 

o„, Br, (r' + ^a") 



^ 2n(2n + l) A„ g„ 

a„+i ~ B„+i (r' + ^^^cr") (r' + na") ~ A2„A2„+i n B„+i 



that after introducing in (C3) we get a simplified version 



(C3) (^{s)y'n (s) = -— ^r"*" "^ Tn (s) y„ (s) - (s) 



Finally using the recurrence relations (C2) we get 

(C4) = 



A„ 2n + 1 A2„ , , , 

« A2n+1 in 



yn (s) + (■5) 



Formulas (C3) and (C4) can be used to calculate solutions of the differential 
equations, fn fact, we put n = in (C4) and we get a differential equation to 
calculate yo{s). 

Taking this value in (C3) we obtain by iteration all the polynomials satisfying 
(CI). 

We can implement these formulas in the discrete case. A polynomial of hyperge- 
ometric type of discrete variable satisfies the following fundamental equations. 
1. Difference equation: 



(Dl) 



a {x) AVy„ {x) + T {x) Aj/n {x) + XnVn {x) = 



where a [x) and t (x) are polynomial of at most second and first degree, respetively, 
and the forward (backward) difference operators are: 



A/(x) = /(x + l)-/(s) , Vf{x)=f{x)-f{x-l) 
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This difference equation can be written in the form of an eigenvalue equation 

A [a {x) p {x) Vyn (x)] + \p {x) yn {x) = 
where p {x) is the weight function satisfying 

A [cr {x) p (x)] =t{x)p (x) 

and 

An = -uAt [xj A a [xj = -n [ T H — a 



is the eigenvalue corresponding to the function y„ (a;) . 

2. Orthogonality relations: 
The polynomial P„ (x) of hypergeometric type satisfy the following orthogonal rela- 
tions: 



b-l 



^ P„ (x) Pm (x) p (x) = dl6r, 



2e = a 



where 5mn is the Kronecker symbol and d„ is some normalization constant. 
3. Three term recurrence relation: 

(D2) xPn (X) = anPn+1 (x) + /3„P„ (x) + In^n-l (x) 



with a„, j3n, 7ra some constants. 

4. Raising and lowering operators: 



(D3) a(x)VP„(x)= ^" 



nr. 



Tn (s) Pn (x) - ^"^Pn+l (x) 



where 



T^M (*■) = T (x + n) + (7 (x + n) — (7 (x) 

Ar„ (s) = At (x) + nA^a (x) 
or T „ = T +na {x) = —- 



2n+l 

because a (x) and t (.x) arc functions of at most second and first degree, respectively. 
We can modify formula (D3) as we did in the continuous case with the help of the 
identities 

n-l . ^ s 

an = Bn n W + -{n + k-l)(j"\,ao = Bo 



an 



fc = 

a„ 2n (2n + 1) n B„ 



Introducing these identities in (D3) we get a more simplified version: 

(D3) a (x) VPn (x) = (^) (^) _ ^a^Pn+i (x) 

n A2„+i 2n 
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This expression defines the raising operator Pn+i {x) in terms of the backward differ- 
ence operator. 

From this expression we can derive another lowering operator in terms of the forward 
difference operator. We substitute the operator V in (D3) by its equivalent difference 
operator V = A — AV and then we introduce the difference equation (Dl) and the 
recurrence relation (D2) obtaining 



(D4) {a{x)+T{x))APn{x) = 



n A2ra+i 2n 2n 

A2n „ 

— — 7„P„_i 
2n 



Pn {X) + 



The advantage of expressions (D3) and (D4) is that all the coefficients are tabu- 
lated. 

As in the continuous case from (D4) , putting n = we get Pq [x) , and inserting 
this value in (D3) we obtain the solutions of (Dl). 

Prom the orthogonal polynomials of hypergeometric type we can construct the 
corresponding normalized functions (up to a phase factor) 



1p„ (x) = rf„ ^ ^/p{x)Pn {x) 

that satisfy equivalent relations and we denote them (NDl), (ND2), (ND3), (ND4) 
in the discrete case and (NCI), (NC2), (NC3), (NC4) in the continuous case. For 
instance, (NCI) becomes: 

a(s)V;'(s) + <7'(s)V;(s) - p(s)-5 (a{s) (p(s)^)') Vn(s) + XnMs) = 



which corresponds to a self-adjoint operator of Sturm-Liouville type. 

3. The Hermite and Kravchuk polynomials. We apply the results of section 
2 to some orthogonal polynomials of continuous and discrete variable and to the 
corresponding normalized functions. 

For the Hermite polynomials Hn (s) we have 

(CI) H"„ (s) - 2sH'^ (s) + 2nH (s) = 

(C2) sHn is) = iif„+i (s) + nff„_i (s) 

(C3) Hn+i (s) = 2sH„ (s) - H'n (s) 

(C4) lf„_i (s) = i^7J'„ (s) 

Introducing the orthonormalized functions 

V„ [s) = {2^n\V^)-^'\-'"I^Hn (s) 

we get 

(NCI) V"„ (s) + (2n + 1 - s^) V'n (s) = 

(NC2) 2sV„ (s) = V2 (n + l)^„+i (s) + v^Vn-i (s) 

(NC3) \/m^n+i (s) = 7f - (^) 

(NC4) V^^n-i (s) = ^ + ^) ^» (^) 
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(NCI) describes the quantum harmonic oscihator, (NC3) and (NC4) are the re- 
aUzation of the famihar creation and annihilation operator 




From (NC4) with n = wc obtain ipQ (s) and inserting this value in (NC3) we 
obtain by iteration the solutions of the harmonic oscillator 

V'n(s) = ^(a"^)"V'o(s) 
V n! 

In the discrete case we take the Kravchuk polynomials kn^ {x, N) with a; = 0, 1, 2, • • • 
N-1 

a (x) = X 

T {x) = {Np — x) /q Tn {x) = {Np — X — n) /q +n 
An = n/q 

an = n+l/q l3n = n + p{N - 2n) , 7„ = (AT - n + 1) 

Inserting these values in the fundamental formulas we get 

(Dl) p{N-x)kn{x + l) + \p{n + x- N)+q{n-x)] kn (x) + qxkn (a; - 1) = 

(D2) xk„ (x) = {n+l) k„+i {x)+[n + p{N- 2n)] kn {x)+pq {N-n + 1) k„-i {x) = 

(D3) (n + 1) kn+i {x) =p{x + n- N) k„ {x) + qxkn {x - 1) 

(D4) q{N-n + l) {x) = {x + n- N) kn (x) + {N - x) fc„ {x + 1) 

For the normalized functions we take the Wigner functions, that appear in the 
representation of the rotation group, d^^^, 

i-l)"'di^, (/?) = d-'^/^)k(r^ {x,N) 

with N = 2j m = j — n m' = j — x p = sin^ ^ ^ q = cos^ ^ 
After substitution we get 

(NDl) VPQiN-x) {x + l)d^_„^ {(3) + 

+ {p{N-x-n) + q{x-n)) 4_„_ + 

+ ./pqx{N-x + l)d^_„^ iP) = 

(ND2) [-p(N -x~n)-q{n- x)] 4„„_ ^_,+ 

+ (n + 1) {N-n)d]_^_^^ {(3) + 

+ ^pqn {N-n+ l)cij_„+i, iP) = 
(ND3) ^pqin+l){N-n)d^_^_,^ = 

= piN-x-n) {(3) + ^pqx{N-x + l)d^_^^ 

(ND4) Vpqn{N-n+l)d^_„^,^ (/3) = 

= p{N-x-n) + ^/pq{x + l) {N - x)d^_^^ (/?) 
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The last four equations can be written down in terms of the new parameters 
j = N/2, m = j — n, m'=j — x, p = sin^/J/2 , g = cos^/3/2 , y^=|sin/3 

(NDl) V(i + m')(j-TO' + l)Cr„'-i + 



+ VU-m') (j + m' + lK™_„,+i (/?) = 

(ND2) V(i + m)(i-m+l)rfi_,^„, (/3) - 

2 



sin /3 



K-m cos (/?) + 



+ - m) (j + m + (/3) = 

L 

2 



(ND3) I sin/3V(i+m)(i-m + lK„_i_„, (/?) 



:sin2^(m + m')d^,„, + 
1 



+ 2 - m') {j + m' + 

(ND4) 1 sin /3v/(i-m)(i + m + (/?) = 



sin2|(m + m')d^,„, + 



+ 2 sill /3 V(J + m') (j - m' + 

Note that (NDl) and (ND2) arc equivalent if wc interchange m m' and take 
in account the general property of Wigner functions 

4,™'(/3) = (-l)"'4',™(/3) 
The same property of duality applies to (ND3) and (ND4). 

In Reference [1] we have constructed creation and annihilation operators with the 
help of (ND3) and (ND4) 



= V(J +m){j-m+ 1) , A =-j= -m){j + m+ 1) 

that together with = ^m, when apply to the spherical functions satisfy the 
50(3) algebra 

[A+,A-]=2A'> , [A^,A'>]=±A^. 

From the connection between the Wigner functions and the solutions of the quan- 
tum harmonic oscillator we have proved [1] the limit relations 

(ND3) ^ (NC3) , (ND4) ^ (NC4) 

Similar results have been obtained by Bijker et al. [6] for the connection between 
the su{2) algebra and the one dimensional anharmonic (Morse) oscillator, and by 
Atakishiev [3] for the lattice implementation of the quantum harmonic oscillator. 
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3.1. The Laguerre and Meixner polynomials. Using tlie general formulas 
of section 2, we get for the Laguerre polynomials L" (s) of continuous variable 

(CI) sL"n (s) + (1 + a - s) L'n (s) + nL„ (s) = 

(C2) (n + 1) i„+i (s) + (n + s) L„_i (s) = (2n + a + 1 - s) L„ (s) 

(C3) sL'„ (s) = - (1 + a - s) i„ (s) + (n + 1) L„+i (s) 

(C4) sL'„ (s) = 2nL„ (s) - (n + a) (s) 

For the normalized Laguerre functions 



we obtain the following differential equation, recurrence relations and expressions for 

the raising and lowering operators: 

(NCI) sV" (s) + s-f/-; (s) + i [-s2 -a^ ^^as^ s\ V^„ (s) + sn^ri («) 



(NC2) v/(n+ 1) (n + a + l)i/'„+i (s) + (n + a)i/'„-i (s) = (2n + a + 1 - s) ^ («) 



(NC3) v/(n + 1) (n + a + l)V«+i («) = ^ (2n + a + 2 - s) V„ (s) + s^; (s) 



(NC4) (n + a)^/.„_i (s) = i (2n + a - s) - (s) 

The first equation, divided by s, corresponds to the self-adjoint operator of the 
Sturm-Liouvillc problem for the normalized Laguerre functions. 

The last two equations can be considered the creation and annihilation operators 
for the normalized Laguerre function. In fact, from (NC4) with n=0 we obtain ■(/'q . 
And from (NC3) we easily obtain 

^n!(Q; + l)„ 

where 

^+V^«(s) = ^/(n + \){n + a + X)■>\)n^x{s) 



A tpn{s) = n{n + a)il)n-i{s) 

The Laguerre creation operator was given by Szafraniec [8] . His formula is equivalent 
to ours if we substitue (NCI) into (NC3). 

Now we take the Meixner polynomials (a;) of discrete variable x, and apply 
the general expressions of section 2: 

(Dl) /X (x + 7) run (x + 1) + (a; - 1) m„ {x - I) + [-fi (a; + n + 7) + n - x] m„ (a;) = 
(D2) urun+i (x) + n (n + 7 - 1) m„_i {x) = [x - I) + n -\- ji {n + 7)] m„ (x) 
(D3) tirun+i (x) = /X (x + n + 7) to„ (x) - xm„ (x - 1) 
(D4) n (n + 7 - 1) m„-i (x) = /li (x + n + 7) m„ (x) - (7 + x) m„ (x + 1) 
For the normalized Meixner polynomials 



MCr) (X) - M"(i-M); M-r(7 + x) . . 

^"^^-V n!(7)„ Vr(^ + l)r(7) ^ 
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we have the following difference, recurrence equations and raising/lowering opera- 
tors: 

(NDl) il + x) {x + l)Mn (x + 1) + ^/^^x{x + -f-l)Mn (x - 1) - 

- [pi{x + n + j) -n + x]Mn{x) = 
(ND2) (7 + n) (n + l)M„+i {x) + (n + 7 — l)M„_i (x) — 

- [a* (x + n + 7) - X + n] M„ (x) = 
(ND3) (7 + n) (n + l)M„+i (x) = ^ (x + n + 7) M„ (x) - 

- ^//zx (x + 7 - 1)M„ (x - 1) 
(ND4) ^ijn (n + 7 - l)M„_i (x) = (x + n + 7) M„ (x) — 

- ^/JlJ^iTx)JxTT)M„ (x + 1) 

As in the case of Wigner functions there exists for the Meixner normalized functions 
an operation of duality between (Dl) <-> (D2) and (D3) (D4) after interchanging 

X <-> n. 

The last two equations can be consider the creation and annihilation operators 
for the normalized Meixner functions: 

(ND3) A+M„ (x) = VM(7 + n) (n + l)M„+i (x) 



(ND4) A- Mr, (x) = (n + 7- 1)M„_i (x) 

From A' Mo (x) = and (ND4) we obtain Mq (x) and taking this value into (ND3) 
we get 

Mn{x)= . (A+)"Mo(x) 

The operators A+, A~ are adjoint conjugate each other and, together with 

A° =/i(2n + 7) 

they close the su(l, 1) algebra 

[A+,A-]^A'> , [A±,A°]=±2M^ 
Finally, from the limit relation between the Meixner and Laguerre polynomials 

it can be proved the limit between the corresponding normalized functions 

and the limit between the recurrence relations 

(ND2) — . (NC2) 

The Meixner polynomials of discrete variable can be used in the solutions of the three 
dimensional harmonic oscillator in conection with the energy eigenvalues [Ref. 7, 
formula (5.10)]. 



CREATION AND ANNIHILATION OPERATORS.. 



9 



The Laguerre-Meixner creation and annihilation operators were presented by F.H. 

Szafraniec at Workshop on Orthogonal Polynomials in Mathematical Physics, June 
24-26, 1996, Universidad Carlos III de Madrid, Leganes, Spain, (but not published). 

Atakishiyev has given recently [10] the hamiltonian and the creation and anni- 
hilation operators for the Meixner oscillator. In fact, his formula (43) and (46) are 
equivalent to ours (NDl), (ND3), (ND4) for the normalized Meixner functions. They 
satisfy the same su{l, 1) algebra if we identify A+ K+, A~ <-> K-, —Kq 

The Laguorrc polynomials of continuous variable arc used in the solutions of the 
hydrogen atom, with different weight function in the orthogonality relations [Ref. 7, 
formula (5.10)]. 

The wave equation for the hydrogen atom can be considered the differential equa- 
tion of the Sturm-Liouville problem for the Lagucrrc functions 

Ms) = €'^Hs) = ^W)L^:^lM 

where 

pi(s)=a(s)p(s)=s2'+2e-« 
In fact, after substituion of V'n(s) in (CI) for the Laguerre polynomials we get 

<(«) - [Pl(s)]"' {[Px{s)^)%n{8) + {n-l- l)s-Vn(s) = 

that corresponds to a self-adjoint operator, from which orthogonality relations for 
V'n(s) with weight function can be obtained. 
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